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ABSTRACT 

The impact of steady increase in surface air temperature on climate change is 
a serious topic of today’s discussion. This study deals with the trend analysis 
and time Series analysis and its future forecasts of surface air temperature of 
Kolkata district of West Bengal state, India during 1901-2002. The mean of 
annual average, maximum and minimum surface air temperature (°C) of last 
102 years are respectively 26.73, 31.26 and 22.23 with 1.27 %, 1.28 % and 
1.73 % coefficients of variation. A non-parametric trend test namely the Mann- 
Kendall (MK) Trend Test along with the Sen’s Slope estimator has been used to 
determine the monotonic trend of temperature and the magnitude of trend 
respectively. On the other hand, the Innovative Trend Test is also being 
performed for pattern determination. Both of the Mann-Kendall (MK) Trend 
Test and Innovative Trend Test show that the trend of average, minimum and 
maximum surface air temperatures of Kolkata are increasing (i.e. upwards) in 
nature. The time series forecasts of surface air temperature are done by 
Autoregressive Integrated Moving Average (ARIMA) model. This upward 
trend of surface air temperature has an impact on human life and urbanization 
as well. The impact of climate change and related socio-economic 
development will affect the ability of a nation to achieve sustainable 
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INTRODUCTION 

Climate change is a burning issue of recent days. The earth’s 
climate is continuously changing its behavior over the past 
century due to variation of surface air temperature, 
precipitation and other climatic conditions. The changing of 
surface air temperature pattern is one of the major impacts 
of climate change. Changes in behavior of surface air 
temperature influences the hydrological cycle too. Moreover, 
the Climatic conditions are responsible for the economic 
condition, livelihood as well as_ overall sustainable 
development of a region. Many researchers have studied on 
different themes of climatic conditions. Pattanaik et.al 
(2007) [1] investigated the analysis of rainfall over different 
homogeneous regions of India. The study was done in 
relation to variability in westward movement frequency of 
monsoon depressions. Mondal et.al. (2012) [2] showed the 
trend analysis of rainfall of north-eastern part of Cuttack 
district of Odisha, India. Anie et.al.(2018) [3]have studied the 
pattern of rainfall in Vamanapuram River Basin of Kerala, 
India by using non-parametric statistical trend analysis. The 
rainfall pattern of the districts of West Bengal has been 
studied by Mukherjee (2017) [4].Girma et. Al (2020) [5] has 
studied on recent trends on climatic variability using 
Innovative Trend analysison upper Huai River Basin. 


In view of the above, study of surface air temperature trend 
is a significant and interesting area of research today. 
Changing pattern of surface air temperature is not only the 
important global and national climatic scenario, but also an 
important scenario for a district too. 
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The Kolkata district is located on the east bank of the 
Hooghly river, it is the major commercial, cultural and 
educational center of East India. The district has a Tropical 
wet-and-dry climate. Huge population as well as pollution 
are the major concerns in Kolkata. The Suspended 
Particulate Matter (SPM) level is high as compared to other 
major districts and cities of India, leading to regular smog 
and haze. Severe air pollution in the district has caused rise 
in pollution-related respiratory ailments. 


The main objective of the present study is to detect the long- 
term monotonic significant trend of mean annual surface air 
temperature (average, minimum and maximum) of Kolkata 
district during 1901-2002 and its related magnitude by the 
Mann-Kendall (MK) Trend Test Method and the Innovative 
Trend (IT) Test Method. The Box-Jenkins autoregressive 
integrated moving average (ARIMA) method is performed to 
show appropriate time series model and future forecasts 
(2003 to 2050) of the mean annual surface air temperature 
(average) of the district. 


Study Area 

The Capital of West Bengal is Kolkata which is the only cent 
percent urban district of West Bengal, India. Kolkata district 
lies between 22° 37’ and 22° 30’ North Latitude and 88° 23’ 
and 88° 18’ East Longitude. As per Census 2011, the district 
covers total 185 sq. Km area with population density 24306 
person/sq. Km. The total population of Kolkata district is 
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4496694 with sex ratio 908, child sex ratio 933, literacy rate 
86.30%, work participation rate 39.90 %. 


Materials 

The datasets for mean annual surface air temperature 
(average, minimum and maximum) (°C) of Kolkata district of 
West Bengal state, India during 1901-2002 have been 
collected from India Water Portal, The Indian Meteorological 
Department (IMD). The demographic data are collected from 
the Census of India, 2011. The Statistical analysis is done by 
using “R” (Version 3.6.1) Statistical software. 


Methodology 

Statistical Inference has basically two parts, one is statistical 
estimation of parameter(s) and other is statistical testing of 
hypothesis. Parametric testing and Non-parametric are the 
two commonly used testing procedures. A parametric 
Statistical test is a test where the model specifies certain 
conditions about the parameters of the population. In this 
case there must be a parent probability distribution form 
which the samples are to be drawn. But in reality, many 
situations arise where the assumptions about the parent 
probability distribution is known to us. Then we have to 
introduce non-parametric approach. A non-parametric test is 
a test where the model does not specify conditions about the 
parameters of the population or the parent probability 
distribution. But the most interesting fact is, in most of the 
non-parametric model, the test statistic follows a Normal 
Distribution asymptotically. 


The mean annual surface air temperature (average, 
minimum and maximum) (°C) of Kolkata district of West 
Bengal state of India during 1901-2002 have been taken into 
consideration separately for analysis purpose. The time 
series analysis is done only on the mean annual surface air 
temperature (average) of the district during 1901-2002 and 
its future prediction have been forecasted for next 48 years, 
i.e. up to year 2050. 


Mann-Kendall (MK) Trend Test Method 

The Mann-Kendall test is a non-parametric test which is 
widely used to detect monotonic trends in series of 
environmental data, climate data or hydrological data. This 
testis the result of the development of the test first proposed 
by Mann (1945) which is further studied by Kendall (1975). 


In Mann-Kendall test we are to test the null hypothesis (Ho): 
There is no significant trend in surface air temperature 
against the alternative hypothesis (H:): There is a specific 
significant trend in surface air temperature at 5% level of 
significance. 


The Mann-Kendall (MK) test statistic is denoted as S and is 
given by 


n-l on 
S=> > sgn(X, — X;,) (1) 

k=l j=k+n 
where, n is the number of data points; Xx; : XxX, are the data 
values for the time period j and k 
respectively; sen(X Fue .4 _) is signum function and is equal 
to 1, 0 or -1 if(X,—X;,) is positive, zero or negative 
respectively. 


Under the null hypothesis, the mean of S is zero and the 
variance of S is denoted by Var (S)and is given by 


Var(S) = in(n—1(Qn+5)-St jt, -DQt,+5}3/18 


where m is the number of the tied groups in the data set and 


t ; is the number of data points in thejtied group. For 


n = 1Othe statistic S is approximately normally distributed 
with mean 0 and variance Var(S). The normally distributed 
test statistic is defined by 


S=l x 
Var(S) ae (3) 
Z=,50 if S=0 
S+1 


4/Var(S) 


The Kendall’s tau is denoted by7 and is given by 





if S<0 


s (4) 
1 LY v2 l _yq/2 
me D-5 LG —DI Pa 1)] 


T= 


In situations where there is no tie(s) occur, then Var (S) and 
T reduce to 


Var(S) = {n(n—1)(2n+5)t/18 (S) 

and 

Sin ne om 
LB n(n —1)] 

respectively. 


Kendall's tau (t) has been also calculated to show the nature 
of trend. The positive and negative 7 value define the 
increasing and decreasing trends respectively. 


Now, for determining the hydro-meteorological magnitude 
of trend the Sen’s slope estimator is used. 


A set of linear slopes is calculated as follows 
dc 2gS* 


i 
Sar (7) 


for(l<k <j <n), dis the slope and X ,, X, are the data 


values for the time period j and k respectively. 


The Sen’s slope( /’) is calculated as the median from all 


slopes. The percentage change (%) is defined by the 
following formula. 


Percentage change = B x Time Length x 100 / Mean (8) 


The positive B value indicates upward (i.e. increasing) trend 
and the negative £8 value indicates decreasing (i.e. 
downward) trend. The test also derives the p-values which is 
the probability of obtaining the observed results of a test 
assuming the null hypothesis is true. This p-value and Z 
value play the significant role for taking decision about 
accepting or rejecting the null hypothesis. 


Innovative Trend (IT) Test Method 

Innovative Trend (IT) Testis widely used technique for 
detection of monotonic pattern of hydro meteorological data. 
In Innovative Trend method, the datasetis classified into two 
classes and then the data points arranged independently in 
increasing order. 
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Let Xi, i=1,2, 3,....n/2 and Xj; j= (n/2)+1,...,.n be placed on X- 
axis and Y-axis on a coordinate system respectively. . If the 
time series data on a scattered plot are collected on the 1:1 
(45°) straight line, it indicates no trend. However, the trend 
is increasing when data points accumulate above the 1:1 
straight line (i.e. upper triangle) and decreasing trend when 
data points accumulate below the 1:1 straight line (i.e. lower 
triangle). 


The trend indicator (@) is defined by 


| 10 (xX ,- X ,) 


where n is the total number of data points, yp is the mean of 
the data series in the first half subseries class. 


i (9) 


The positive @ values defines the increasing trend and 
negative ~ value defines the decreasing trend. When the 
scattered data points closest around 1: 1 straight line, it 
indicates non-existence of trend. 


Box-Jenkins Autoregressive Integrated Moving Average 
(ARIMA) Method: 


A time series is defined as an arrangement of data over 
specific time period [6]. Box-Jenkins autoregressive 
integrated moving average (ARIMA) method due to Box and 
Jenkins. ARIMA is a very popular method used to produce 
accurate future forecasts. This ARIMA (p, d, q) method is a 
combination of autoregressive (AR) of order p and moving 
average (MA) of order g with degree of differencing involved 
(d). The differencing (d) is used to make the data stationary 
to non-stationary. Therefore, the model is called as 
Integrated ARMA, i.e. ARIMA. 


Let us assume Y;, t=1(1) n bea time series then the general 
form of difference taking time series is given by 


d d-l d-l 
VY =VY,-V"Y_, (10) 
where d is the degree of differencing involved. 


An autoregressive moving average (ARMA) model with 
order of autoregressive part (p) and order of moving average 
part (q) is defined by ARMA (p,q) model and the model is 
expressed by the following model expression 


Results and Discussion 


p Y 
y, =u+) OY, ; +) PE, TE; (11) 
i=l i=l 


where J/ isthe mean value and 9&@ are the parameters of 
autoregressive and moving average part respectively which 
are to be estimated. Itis to be noted that the errorterm(é€, ) 


is the white noise which is i.i.d (0, o ). A time series is said 
to ARIMA (p, d, q) if V“Y, is ARMA (p, q). 


The Augmented Dickey Fuller (ADF) test is used to check if 
the time series is stationary or not. The p-value of the ADF 
test defines whether the time series is stationary or not 
where the null hypothesis is taken as the time series is not 
Stationary against the alternative hypothesis is that the time 
series is stationary at a pre-assigned level of significance. For 
making the time series stationary, the method of differencing 
is needed. 


Autocorrelation function (ACF) graph and_ partial 
autocorrelation (PACF) graph are used to find the initial 
degree of qand pin ARIMA (p,q) model. ARIMA model deals 
with stationary data. 


Based on the values of p, d and gq, some special cases are 
discussed as follow: 


White noise model: ARIMA (0, 0, 0) with p=0, d=0, q=0 


Random walk: ARIMA (0, 1, 0) with no constant with p=0, 
d=1, q=0 

Random walk with drift: ARIMA (0, 1, 0) with constant with 
p=0, d=1, q=0 

AR (p): ARIMA (p, 0, 0) with d=0, q=0 

MA (q): ARIMA (0, 0, q) with p=0, d=0 

ARMA (p, q): ARIMA (p, 0, q) with d=0 

Akaike’s information criteria (AIC) isapplied for selecting the 
appropriate ARIMA model. Usually the model with lowest 


value of AIC among the all possible models is to be selected 
as the best fitted model. 


Based on the best fitted ARIMA (p, d, q) model, forecasts of 
the time series isobtained with 80 % and 95 % (low and 
high) prediction intervals. 


The mean annual surface air temperature (average, minimum and maximum) during 1901-2002 have been taken into 


consideration. 


Table 1: Descriptive Statistics Table for Mean Annual Surface Air Temperature of Kolkata District 


Temperature 


27.76 26.73 


Minimum Maximum Mean Median 


Temperature 


26.73 


Coefficient of variation (%) Skewness_ Kurtosis 


1.27 - 0.11 


Temperature 


Temperature 





N.B.: The table is tabulated by the author 


Table 1 shows that the mean of annual average, maximum and minimum surface air temperature (°C) of last 102 years are 
respectively 26.73, 31.26 and 22.23 with 1.27 %, 1.28 % and 1.73 % coefficients of variation and median values 26.73, 31.24 
and 22.18. The skewness values define that the annual average (-0.11) and maximum (-0.08) temperature are skewed left in 
nature where minimum temperature is skewed right (0.33). The kurtosis values for annual average (0.69) and maximum (0.87) 
temperature are positive and the for the case of minimum ( -0.15) temperature the kurtosis value is negative during 1901- 


2002. 
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Table 2: Trend Test Table for Mean Annual Surface Air Temperature of Kolkata District 
Mann-Kendall Trend Test Innovative Trend Test Significant 


ee Kendall’s Sen's Se Innovative pat 
Temperature Z Kendall % Trend (at 95% 
ae Tau Slope p-value eee Trend 
value sae (1) (B) Change ice Confidence 


ope) Interval) 





Temperature 
Temperature 
Temperature 


N.B.: The table is tabulated by the author 


From Table 2, it is clear that for mean annual surface air temperature (average, maximum and minimum), the Z-values,Mann- 
Kendall Statistic (S) and Kendall's tau (t)values are all positive. The Sen’s slop estimator (() values are also positive for all the 
above three cases. The MK testis performed at 5 % level of significance. The p-values are all less than 0.05. Hence it is clear that 
the trend of meanannual surface air temperature (average, maximum and minimum) are all significantly increasing (i.e. 
upward) at 5% level of significance or 95% confidence interval with percentage change 1.95 %, 1.37% and 2.57% during 1901- 
2002. In case of Innovative Trend test, the innovative trend indictors (q@) are all positive and hence it may be found that the 
trend for all the above three cases are significantly increasing. From this it may be found that for both trend test viz. Mann- 
Kendall Test and Innovative Test, the mean annual surface air temperature (average, maximum and minimum) of Kolkata 
district are all significantly increasing in nature during 1901-2002. 


A time series is defined as an arrangement of data over specific time period. ARIMA model is widely used reputed time series 
model for forecasts or future prediction. 


The mean annual surface air temperature (average) during 1901-2002 have been taken into consideration for time series 
analysis. The data is transformed to a time series data to carry out the ARIMA model analysis. Then we are to check whether the 
data is stationary or not. The stationary is checked by the Augmented Dickey Fuller (ADF) test at 5% level of significance. The 
p-value plays a significant role for taking decision about accepting or rejecting null hypothesis. The test result is tabulated as 
follows: 


Table 3: Augmented Dickey Fuller (ADF) test 
Test Order of differencing (d) involved p-value Remarks 


PADF/] dO | 0.64 _| Non Stationary | 





N.B.: The table is tabulated by the author 


From Table 3, itis observed that the 1° order differencing of the time series is stationary with p-value 0.01 (< 0.05) as the test 
is performed at 5 % level of significance. 


Now by using “R” Statistical software (Version 3.6.1), the ARIMA model is constructed for 1st order differenced (d=1) data 
series. 


Table 4: ARIMA Model Summary Table 
Model u AR(1) AR(2) AR(3) AR(4) MA(1) MA(2) MA(3)_ AIC 





| ARIMA (4,1,3) | 0 | -0.6560 | -1.2055 | -0.1343 | -0.2979 | 0.0184 | 0.5385 | -0.6702 | 43.49_ 


N.B.: The table is tabulated by the author 


From Table 4, it is found that for ARIMA (4,1,3) with zero mean, the calculated value of Akaike’s information criteria (AIC) 
(43.29) is the smallest among the possible models. Hence the model “ARIMA (4,1,3) with zero mean’ is selected as the best 
fitted ARIMA model for mean annual surface air temperature (average) of Kolkata district during 1901-2002. 


The model expression of ARIMA (4,1,3) with zero mean is given by the following equation: 


Y! =-0.6560Y,, -1.2055Y!, —0.1343Y,, —0.2979Y!, +0.0184e,_, + 0.5385¢,_, -0.6702¢,, +, ao 


where Y/, is the 15t order difference of the time series Y 


Based on the above ARIMA (4,1,3) with zero mean, the future forecasts for mean annual surface air temperature (average) of 
Kolkata district from 2003 to 2050 is done with 80% and 95% (low and high) prediction intervals. 
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Figure 1: Forecasts from ARIMA (4, 1, 3) 


Figure 1 depicts the actual time series data on mean annual surface air temperature (average) of Kolkata district during 1901 
to 2002 and its forecasts from ARIMA from2003 to 2050 at 80% and 95% (low and high) prediction intervals. 


Concluding Remarks 

The impact of climate change on annual surface air 
temperature is a serious issue now a days. The mean annual 
surface air temperature (average, minimum and maximum) 
are all significantly increasing in nature. The time series 
forecasts from ARIMA of mean annual surface air 
temperature (average) is not also decreasing in pattern. For 
a district with huge population and pollution like Kolkata, 
the increasing trend of surface air temperature is a very 
alarming agenda today. Rapid growth of population, 
concretization, urbanization, factories & industries, emission 
of greenhouse gases, destroying trees and green places and 
other alike activities are mainly responsible for long-term 
upward trend of surface air temperature of Kolkata. 
Monotonic upward pattern of surface air temperature has a 
deep impact on environment, livelihood, health, socio- 
economic conditions and overall sustainable development of 
a region. This study will give the overall information to the 
Government, community, planners, researchers, policy 
makers about the long-term increasing pattern of surface air 
temperature of the district and help them to make necessary 
action plan and framework to save our environment and 
planet from the climatic disasters in near future. 
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